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Abstract 

The properties of spin polarized neutron matter are studied both at zero and finite temperature 
using Skyrme-type interactions. It is shown that the critical density at which ferromagnetism 
takes place decreases with temperature. This unexpected behaviour is associated to an anomalous 
behaviour of the entropy which becomes larger for the polarized phase than for the unpolarized 
one above a certain critical density. This fact is a consequence of the dependence of the entropy 
on the effective mass of the neutrons with different third spin component. A new constraint on the 
parameters of the effective Skyrme force is derived if this behaviour is to be avoided. 

PACS numbers: 21.30.-x, 21.65.+f, 26.60.+C, 97.60.Jd 
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I. INTRODUCTION 

Since the pioneering work of Vautherin and Brink Q there has been an intensive use 
of Skyrme effective nucleon-nucleon (NN) interactions to study properties of finite nuclei 
and nuclear matter, the latter one mainly in conditions of astrophysical interest. This type 
of phenomenological NN interaction is thought to be used within a Hartree-Fock scheme 
(HF). Its zero-range character leads to simple analytical expressions for basic properties of 
symmetric nuclear matter such as the binding energy (a v ), the saturation density (po), the 
incompressibility modulus (K^) or the symmetry energy (a s ). Actually, the experimental 
values of these quantities together with the binding energy of some doubly magic nuclei have 
been traditionally used to fit the parameters entering the general expression of a Skyrme 
force. The main advantage of these forces comes from their analytical character, which 
makes them very useful to get a physical insight into problems where the fully microscopic 
calculations are either very time consuming or not yet possible to implement. Once a 
parametrization is determined, one can, for instance, study nuclei and nuclear matter under 
conditions far from those used to fix the force in quite a simple way. In doing so, however, 
one must always keep in mind the particular limitations of the Skyrme parametrization 
which is being used. 

By construction, most of the Skyrme forces used in the literature are well behaved around 
the saturation density of nuclear matter and for moderate isospin asymmetries. However, 
not all the Skyrme parameters are completely well determined through the fits of given sets 
of data and only certain combinations related to the basic properties mentioned above are 
really empirically determined This leads to a scenario where, for instance, different 
Skyrme forces produce similar equations of state for symmetric nuclear matter but very 
different results for neutron matter. This is easily understood if one considers that neutron 
matter or, equivalently, the systems with large isospin asymmetries, are not part of the 
common input data that determine the parameters of the interaction. Obviously, this feature 
should be corrected if Skyrme-type forces are to be used in conditions of large neutron 
to proton ratios such as nuclear matter inside neutron stars or nuclei near the drip line. 
Recently, several sets of Skyrme forces have been constructed taking into account different 
data coming from highly isospin asymmetric systems. The most well-known among these 
are probably the parametrizations of the Lyon group (SLy interactions) 0,0, which also 
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iook into account variational results for neutron matter obtained with a realistic interaction 
ij] to fix the parameters of the force. Another important piece of input for neutron-rich 
systems are the isotope shifts of medium and heavy nuclei. In that case, modifications of 
the spin orbit term of the force have been taken into account to correctly reproduce the 
data. With these prescriptions, the Ski parametrizations were created almost a decade ago 
[f|. Recently, an extensive and systematic study has tested the capabilites of almost ninety 
existing Skyrme parametrizations to provide good neutron-star properties [?|. It was found 
that only twenty seven of these forces passed the restrictive tests imposed, the key property 
being the behaviour (increasing) of the symmetry energy a s with density. 

Another situation of astrophysical interest refers to the possibility of the spontaneous 
appearance of spin polarized states in nuclear matter. This type of instability, i.e., a fer- 
romagnetic transition at high densities, has been studied since long ago using different 
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nevertheless, still contradictory. On the one hand, Skyrme interactions predict several types 
of instabilities at increasing density Q]. In particular, currently used Skyrme forces show 
a ferromagnetic transition for neutron matter at densities in the range (l.l-3.5)po (22I 23^ . 

protons into the 
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However, it has recently been shown that by including a small fraction o 
system, the onset density of ferromagnetism can be substantially reduced 
land, recent Monte Carlo simulations j^j and also Brueckner-Hartree-Fock calculations 
2fiL l27| using modern two- and three-body realistic interactions exclude such an instabil- 
ity, at least at densities up to five or six times po- This transition could have important 
consequences for the evolution of a protoneutron star, in particular for the spin correlations 
in the medium which do strongly affect the neutrino cross sections and the neutrino mean 
free paths inside the star Therefore, drastically different scenarios for the evolution of 
protoneutron stars are to be considered if the existence of such a ferromagnetic transition is 
confirmed. 

Most of the studies of the ferromagnetic instability have been conducted at zero tempera- 
ture. However, the description of protoneutron stars requires a study at finite temperature. 
Thus, since the general conditions for ferromagnetism are well established in the case of 
the Skyrme forces at zero temperature and even though their predictions are quite different 
than those of microscopic calculations, it is still interesting to study how the situation is 
modified at finite temperature. Intuitively, due to the thermal disorder, one would expect 
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the onset density for ferromagnetism to increase with temperature. But, as we will see later 
on, the behaviour turns out to be the opposite one and what actually happens is that the 
ferromagnetic transition takes place at smaller densities. In this work, we associate this fact 
to an anomalous behaviour of the entropy as a function of the spin polarization. However, it 
is worth noticing here that there is no thermodynamical inconsistency in this dependence. If 
eventually, microscopic calculations with realistic interactions would not confirm this trend, 
we find out a new constraint on the parameters of the effective Skyrme force in order to 
avoid such a behaviour. 

In the next section we summarize the expressions of the energy and the single-particle 
energies at zero temperature and find the free energy and entropy at finite temperature for 
neutron matter as a function of the spin polarization. In section III, we present some results 
involving the thermodynamical potentials of the system. Section IV is devoted to the study 
of the entropy, its relation to the effective mass and the possible constraints which could 
be imposed to the parameters in order to avoid such an anomalous behaviour in both the 
classical and the low temperature limit. Finally, the main conclusions are summarized in 
section V. 

II. POLARIZED NEUTRON MATTER 

Most of the Skyrme interactions used in the literature have the following general form: 
V(r 1 ,r 2 ) = to(l+x P a )S(r) + lt 3 (l + x 3 P a )[p(R)} a 5(r) 



with r = r x — r 2 , R = (r x + r 2 )/2, k = (Vi — V 2 )/2i the relative momentum acting on 
the right and k' its conjugate acting on the left. P a — (1 + U\ ■ cx 2 )/2 is the spin exchange 
operator. The last term, proportional to Wo, corresponds to the zero-range spin-orbit term, 
which does not contribute in homogenous systems and thus will be ignored for the rest of 
the paper. 

Let us consider a homogeneous system of neutrons characterized by a total density p, 
which is the sum of the spin- up (pj) and spin-down (pj densities. At zero temperature, the 
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total energy will be a function of pj and pi or, alternatively, of the total density p and the 
spin polarization parameter A, defined as A = (pf — pj.)/p. 

The total energy of this system in the Hartree-Fock approximation is given by the sum 
of the kinetic energy associated to the Fermi gas of polarized neutron matter and the expec- 
tation value of the Skyrme interaction between the wave function describing two free Fermi 
seas corresponding to neutrons with two different spin orientations. These two Fermi seas 
have all the single-particle states characterized by a good linear momentum, occupied up to 
the Fermi levels k F ^ and kpi defined through kpm) = (6n 2 p^^y^ 3 . 

The total energy per particle for zero temperature has the following expression: 

h 2 l l 

e (PhPl) = b~t + T l\ + J- p [ 2t 2(l + x 2 )\ [r T p T + T iPi \ 

+^ [*i(l - x x ) + t 2 (l + x 2 )\ [r lPl + r lPl ] 

IT/ x 1 / 
+- t (l-x ) + -t 3 (l-x 3 )p 



P\Pi ■ (2) 



P 

To simplify the notation, we will use from now on the symbol a to indicate the third spin 
component. The functions r CT are related to the average kinetic energy of the Fermi model 
of polarized neutron matter: 

r, = jj(67rV) 2/ V = ^ (3vr 2 p) 2/3 p(l ± A) 5 / 3 , (3) 

where the plus (minus) sign corresponds to the up (down) spin projection. From the energy 
per particle one can easily derive the chemical potentials (up and down) and the pressure. 
Another important quantity for our analysis is the single-particle energy: 

h 2 k 2 

ta(k) = —— + U r7 (k,p hPl ) . (4) 

The single-particle potential U a (k, pj , pj takes into account the interaction of a particle with 
momentum k and spin projection a with all the rest. It has a quadratic dependence on the 
momentum which is usually incorporated in the single-particle spectrum as a momentum 
independent effective mass: 

h 2 k 2 

to{k) = — - + U a (p hPl ) , (5) 
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where the effective mass is given by 

rri 
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m l + 2™a a (p hPl )' 
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The momentum independent part U a of the single-particle potential U a is then given by: 

1 



Ua(pr,Pi) 
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It is worth mentioning that the expression for U a contains also the rearrangement term 



Un(p^Pl) = ^at 3 p a 1 p a p-a, 



(9) 



which takes into account the effect of the density dependence of the effective interaction on 
the single-particle potential. Using this prescription for the single-particle potential, one can 
check that the chemical potential calculated from the energy per particle: 



P^Pi) = e(p hPl ) + p{^^ 



(10) 



does exactly coincide with the single-particle energies at the respective Fermi surfaces of 
each species, p a = e(/cF CT ). 

The extension of these expressions to finite temperature is rather straightforward. The 
expression for the internal energy and for the single-particle energy are exactly the same as 
in the zero temperature case, the only change coming from r CT , which at finite temperature 
is given in terms of the so-called Fermi integrals 



g (2m* \ 5/2 , 



(11) 



(2vr) 2 V H 2 , 

where g is the spin degeneracy factor of the system (in our case g = 1 for each spin compo- 



nent) and where: 
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dx- 



/o 1 + e x -v 

The parameter r] a should be calculated by inverting the equation: 
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which simply states, in terms of r\ ai that the integration of the Fermi momentum distribution 
of each spin component: 

n„{k) = (l + e'-^^y 1 (14) 

should coincide with the density of the component. 

Once rj a has been obtained, one can also calculate the chemical potentials at finite tem- 
perature: 

P*(phPh T ) =V*T + U a (p h p h T) (15) 
and the entropy per particle: 

s(phPh T ) = ^j^(p h p h T) + ^ Sl (p h p h T) (16) 

where the entropies of each component are given in the Hartree-Fock approximation by: 

1 f d 3 k 

sAPhPh T ) = / 77T\3 K(A:)lnn CT (A;) + (1 - n CT (A;)) m(l -n a {k))} 

5 1 1 (2m* \ 3/2 

Finally, for a fixed density and temperature, the suitable thermodynamical potential is 
the free energy. We can get if from the previous expressions of the energy and the entropy: 

f(ph Ph T ) = e{phPl> T ) ~ TsiPhPh T ) ■ ( 18 ) 

From the free energy per particle, we can get the rest of the macroscopic properties of the 
system as, for instance, the pressure (and thus the EoS). In our case, we are particularly 
interested in the inverse magnetic susceptibility which can be obtained from a second 
derivative of the free energy with respect to the spin polarization: 

X p?p \dA 2 J A=0 

where p is the magnetic moment of the neutron. 

Notice that in our approach the effective mass m* does not depend on the temperature 
and that the chemical potential obtained from the normalization condition of the density 
of each component (|13|). when the rearrangement is taken into account, coincides with the 
chemical potential derived from the free energy through its derivative with respect to density. 
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III. ENERGETICS OF POLARIZED NEUTRON MATTER 

For the discussion of our results we have chosen the SLy4 ^] and the SkI3 [fj Skyrme 
forces, both of which passed successfully the careful tests of Ref. 0]. In Table HI we report 
the values at saturation density of the binding energy, symmetry energy, incompressibility 
modulus and maximum mass of the neutron star which can be obtained using the EoS of 
these two Skyrme forces when the f3- stability conditions in the presence of electrons and 
muons are imposed. 

Since we are interested in analyzing the behaviour of these forces with respect to the 
polarization of neutron matter, in Fig. (Q) we report the ratio between the inverse magnetic 
susceptibility of interacting neutron matter and that of the corresponding free Fermi gas of 
neutrons as a function of the density for different temperatures. At zero temperature, as it is 
well known, both forces present a magnetic instability, i.e., the previous ratio becomes zero 
at the critical densities (p c = 0.60 fm~ 3 for SLy4, p c = 0.37 fm -3 for SkI3), but, contrary 
to what can be intuitively expected, the onset density for the magnetic instability decreases 
with temperature for both forces. It is precisely this anomalous behaviour that we want 
to explore here. Since we are mainly interested on the study of thermal effects, we have 
explored unrealistically high temperatures (much higher than those needed in the evolution 
of proto-neutron stars) in order to magnify them. 

A complementary information which is rather helpful in this analysis is the difference 
between the free energy of totally polarized and unpolarized neutron matter. This difference 
as a function of density is reported in Fig. (j2J). Once this difference has become negative, the 
totally polarized system will have a lower free energy and therefore the system will prefer 
a polarized phase in front of the non-polarized one. Notice that the density at which this 
difference becomes negative does not coincide with the onset of the magnetic instability. 
The critical densities defined by this criterium at zero temperature are pf = 0.71 fm -3 and 

= 0.44 fm -3 respectively, both of them larger than the corresponding p c 's. The critical 
density p c signals the density at which the unpolarized phase becomes unstable around 
A = 0, which, however, does not imply that the system prefers the fully polarized phase, 
i.e., the minimum of the free energy is not necessarily located at A = 1. Even more, beyond 
pf, as we will check later on, one can not guarantee that the minimum of the free energy is 
located in the fully polarized system. 
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In this discussion, it is worth noticing that the energy per particle for fully polarized 
neutron matter with Skyrme interactions reduces to: 



This can be understood if one considers that the two body states of fully polarized neutron 
matter are all triplet states on both the spin and the isospin spaces, i.e., they are symmetric 
in the spin-isospin variables, and, therefore, due to the Pauli principle, they do not see s- 
wave contributions (associated to the purely contact terms of the Skyrme force), but only 
p-waves originated from the gradient terms. As the parameter £ 2 is usually taken negative, 
one has to take £2 < —1 in order to avoid the collapse of fully polarized neutron matter 
. Both SLy4 and SkI3 fulfill this condition. For SLy4, for instance, = —1. For this 
interaction, totally polarized neutron matter has m*/m = 1 and its energy reduces simply 
to the kinetic energy of a free Fermi sea. 

The free energy is composed by the sum of two contributions, the internal energy and the 
entropy. Therefore, in order to understand the anomalous behaviour of the free energy, it 
is reasonable to analyze how both terms behave as a function of density. With this aim, we 
show in Fig. (j3J) the internal energy difference e(p, A = 1, T) — e(p, A = 0, T) as a function 
of the density for several temperatures. As expected, the density at which this difference 
becomes zero grows with temperature, so the origin of the surprising behaviour detected in 
Fig. (j2j) should be attributed to the entropy contribution of the free energy. 

The difference between the entropy of the polarized and the unpolarized phases is shown 
in Fig. as a function of the density for three different temperatures. One would naively 
expect the entropy of the polarized system to be lower than that of the unpolarized one 
because, intuitively, the fully polarized phase is more "ordered" than the unpolarized one. 
This is, in fact, the behaviour of the free Fermi sea. However, for the interacting system, 
the modification of the single particle properties, in particular those due to the effective 
mass, can invert this behaviour even in the Hartree-Fock approximation. This inversion 
takes place at relatively low densities, pf = 0.15 fm~ 3 and pf = 0.08 fm~ 3 for SLy4 and 
SkI3 respectively, much smaller than the onset densities at which ferromagnetism appears. 
Notice also that the entropy critical densities are temperature independent. 




(20) 



where the effective mass m* depends only on the parameters £2 and %2- 

m*(p, A = 1) 1 



(21) 



m l + f4* 2 (l + x 2 )p ' 
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As we have seen previously, there is a certain range of densities (between p c and pf) 
where neutron matter is unstable around A = but the fully polarized phase is still not 
energetically favourable. In this range (and also for densities a bit higher than p£), one can 
see that the minimum free energy happens at a partial polarization < A < 1. This is 
seen in the left panel of Fig (j5J), where we show the free energy of the system as a function 
of polarization at zero temperature (which of course coincides with the energy per particle) 
and for several densities, obtained with the SkI3 interaction. At densities below p c (the 
critical density provided by the susceptibility criterium), the minimum energy takes place at 
A = (this is the case of the first considered density, p = 0.32 fuT 3 ). When we overpass p c , 
the minimum appears at an intermediate polarization and the energy of the fully polarized 
phase is still larger than the unpolarized one (this would be the case for p = 0.40 fm~ 3 ). 
When the density increases, the minimum moves to larger values of A, and, beyond pf, we 
observe that f(p, A = 0) > /(p, A = 1) (this is the case for p = 0.48 fm~ 3 ) even though the 
A = 1 phase is not the preferred one. Finally, for larger densities, the minimum shows up 
at the fully polarized configuration. 

One can also observe the appearance of a preferred partial polarization as a function of 
the temperature for a given density, as it is illustrated in the right panel of Fig (|SJ) where the 
free energy per particle at p = 0.36 fm~ 3 is drawn as a function of the polarization. At this 
density, the internal energy has the minimum at A = for all the temperatures. However, 
the chosen density is above the onset density for the anomalous behaviour of the entropy 
p > pf, i.e., the entropy is larger for the polarized phase. As a result, the free energy 
develops a minimum at < A < 1 that moves to higher polarizations when the temperature 
increases. Both panels of Fig. (j3J) are typical examples of a spontaneous symmetry breaking, 
where the interaction drives the ground state of the system to non-zero polarizations. 

The spontaneous breaking of the spin rotational symmetry indicates the appearance of a 
phase transition in the system. In our case, this phase transition is between a non-polarized 
and a polarized phase and the natural order parameter is, thus, the polarization. In Fig. 
(JHJ) we report the polarization associated to the minimum of the free energy. For a given 
density and temperature, this is the polarization that gives the lowest free energy and it 
indicates, then, the polarization at which our system is thermodynamically stable. With 
this figure in hand, we can check all the results that we have been discussing up to now. At 
zero temperature, below p c , the equilibrium configuration corresponds to A = 0. Over p c , 
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the polarization grows steeply up to A = 1. When the temperature increases, the curves 
shift to smaller densities, so p c decreases with temperature and also the density at which the 
system becomes fully polarized. Of course, these results are somewhat academic since the 
densities needed for neutron matter to be fully polarized (in particular for SLy4) are very 
high and, thus, unattainable in a proto-neutron star. However, we find them useful since 
they serve to confirm the results concerning the anomalous behaviour of the free energy with 
temperature. 

IV. ENTROPY AND EFFECTIVE MASS 

In this section we investigate the relation between the anomalous behaviour of the entropy 
and the dependence of the effective mass on the polarization. To this end, we will determine 
the density at which the difference between the entropy per particle of the fully polarized 
and the unpolarized phase becomes positive, i.e., does not show the expected behaviour 
s(p, T, A = 1) — s(p, T, A = 0) < 0. We will explore two situations: the classical limit, 
defined by the condition ^ — > (with A = being the de Broglie wavelength) and the 

degenerate limit, where — << 1. 

In the first place, we consider the classical limit. In the interacting case, it is useful to 
introduce the de Broglie wave length associated to the effective mass: 



It is worth reminding that, for Skyrme interactions, the effective mass depends on the density 
but is not affected by the temperature. In addition, it turns out that one can easily write 
all the relevant quantities for the non-polarized system (A = 0) and the fully polarized one 
(A = 1). For instance, the internal energy can be casted in the form: 



where Ui(p,A) is the contribution to the internal energy that can not be included into the 
effective mass terms: 




(22) 



e cla (p,T,A) = ^T + U 1 (p,A) , 



(23) 



U x {p, A) = P - [*o(l - x ) + ^3(1 - x 3 )p a j (1 - A 2 ) . 



(24) 



For the entropy, one has: 



Sda(p,T,g) = - - In 



pX 



(25) 
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where we have explicitly shown the dependence on the degeneracy factor g. In order to get 
the fully polarized case (A = 1), it is enough to set g — 1, while the non-polarized case 
(A = 0) is obtained by setting g = 2. The previous expression can be splitted in two pieces, 
the entropy of a free gas plus a correction term associated to the effective mass: 



s cla (p,T,g) 



5_, n 




3, 
+ - In 


~m*~ 
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. 9 _ 


2 


. m _ 



(26) 



where one has to take into account that the effective mass depends on both the density p 
and the polarization A. For the free case, one can in fact check that the expected inequality 
is fulfilled at all densities and temperatures: 



s cla (p, T, A = 1) — s da (p, T, A = 0) = - In 2 < . 
In the interacting case, however, this difference becomes: 

3 

As da (p) = s da (p, T, A = 1) - Sd a (p, T, A = 0) = - In 2 + - In 



m*(p, A = 1) 



(27) 



(28) 



m*(p, A = 0)_ 

which clearly shows the influence of the interaction, via the effective mass, to the entropy 
difference, and which turns out to be independent of the temperature. If we require the 
difference to be negative, the following condition should be satisfied: 

m*(p, A = T 



< 



2 2/3 _ 



(29) 



m*(p, A = 0) 

Taking into account the expressions for the effective mass of the polarized and unpolarized 
phases for the Skyrme interactions (jOJ, one gets the following condition for the parameters 
of the effective force: 

l + f?| [ti(l-x 1 ) + 3t 2 (l + x 2 ) 



l + f?f2t 2 (l + x 2 ) 



< 



2 2/3 _ 



(30) 



The previous expression is not well defined whenever the denominator becomes zero. In 
fact, one can check that eleven of the twenty seven forces that passed the tests of [t| have a 
vanishing denominator at densities below 0.5 fm -3 . Since this denominator is nothing but 
the inverse effective mass for polarized neutron matter, such a singularity shows that these 
parametrizations were not devised to describe highly spin and isospin asymmetric matter. 
For the rest of the forces, one can check that m*(p, A = l)/m*(p, A = 0) is a monotonically 
increasing function of the density, so we can univocally define a critical density pf a . For SkI3, 
for instance, one can use the condition ()30|) to obtain the critical density p c ^ a = 0.08 fm~ 3 . 
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The case for Skyrme Lyon forces is special, since most of them have X2 = — 1 (the only 
exceptions are SLyO and SLy2). In that simpler case, the previous condition reduces to: 

l + ^|t 1 (l-x 1 )<2 2 / 3 . (31) 

Then, as £i(l — x{) > 0, one can easily obtain a critical density (p^ la = 0.15 fm -3 for SLy4) 
beyond which one can be sure that, if the temperature is high enough to reach the classical 
limit, one will observe an anomalous behaviour of the entropy. In fact, one can check that 
the temperature is not an essential factor since we have seen that pf (precisely the density at 
which s(p,T,A = 1) — s(p, T, A = 0) changes sign) is quite independent of temperature and, 
thus, this classical entropy critical densities coincide with the pf defined before, pf a = pf . 
In Table |H] we give a list of these densities for the sixteen forces which have no singular 
effective mass for neutrons in fully polarized matter and compare it to the corresponding 
ferromagnetism onset densities at zero temperature p c . 

In second place, we investigate the low temperature limit. In this case, we can express the 
different thermodynamic quantities for the polarized and the unpolarized phases in terms 

of the Fermi momenta kp = (67r 2 pj g) 1//3 and the Fermi energies associated to the effective 

h 2 k 2 

mass e F (A) = 2to *(a) • ^ n ^ ne ^ ow temperature regime, the internal energy can be written as 
(this expression is only valid for A = or A = 1): 



3 

ei ow (p,T,A) = ge>(A) 



, 5-7T 2 ( T 

1 + 



12 V^(A) 



+ U 1 (p,A), (32) 



where U\ is the same function defined previously for the classical case and where we recover 
the well-known T 2 dependence of the internal energy. The entropy per particle (again only 
valid for A = or A = 1), on the other hand, shows a linear dependence on T: 



IT 2 _ I 2 



S - ( "- r - A) =2^ T =3/ (0)T - (33) 

where we have introduced the density of states at the Fermi surface, N(0) = grrfkp jli^h 1 '. 

At this point, we can perform the same analysis that we made for the classical limit. Let 
us start again by the free case and rewrite (|3*3*jl to see the explicit dependence on the spin 
degeneracy: 

s low (P, T, A = 1) - s low (P, T, A = 0) = - 2 2 / 3 ) , (34) 
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which in fact is smaller than zero. For the interacting case, the only difference comes, again, 
only from the effective masses: 



Asi ow (p, T) = s iow (p, T, A = 1) - si ow (p, T, A = 0) = 
n 2 Tm*(p, A = 0) \ m*(p,A = l) _ ^ 
_m*(p, A = 0) 



(35) 



h 2 (6n 2 p) 2 / 3 

By requiring this difference to be negative, we recover the same condition for the effective 
masses that we found in the classical regime (jHOJ)- Contrary to the classical limit, in this case 
Asi ow depends on the temperature, although this dependence can not change its sign. Now 
we can understand why is independent of the temperature: since in the classical limit 
(which requires high temperatures) and in the low temperature limit the entropy critical 
densities are exactly the same, one can guess that for intermediate temperature pf will also 
not change. 

We have seen that the interaction influences the entropy only through the effective mass 
and we have been able to trace back the anomalous behaviour of the entropy (as a function 
of the polarization) to the violation of relation ([30)1 . which does only involve effective masses. 
In Fig. (J7j). then, we show the effective mass of neutrons with both orientations of the third 
spin component as a function of polarization for a fixed density p = 0.16 fm -3 , which lies 
above the entropy critical density pf for both SLy4 and SkI3. Obviously, for A = 0, the 
effective masses for the neutrons with spin up and spin down coincide. The horizontal 
line in both panels is placed at m*/m = 2 2 ' 3 m*(A = 0)/m, indicating the upperbound 
for m*(A = l)/m, above which we will find the anomalous behaviour of the entropy. The 
effective mass of the most abundant species increases with the polarization, while the one of 
the less abundant is the decreasing one. Notice that for SLy4, as we have already pointed 
out, the effective mass for the fully polarized system is simply the bare mass. 

The behaviour of the entropy as a function of the polarization at a fixed density p = 
0.32 fm -3 for several temperatures and for both Skyrme forces is shown in Fig. (jHJ). As 
expected, the entropy increases with temperature. However, since this density is above pf, 
the entropy of the fully polarized phase is larger than that of the unpolarized one, giving the 
anomalous behaviour of the entropy with the polarization. A careful look shows that this 
difference increases with temperature as expected from the results of Fig. (jlj). For a given 
temperature and polarization and considering that the entropy in both limits (and also in 
between them) is an increasing function of the effective mass and that, as we have seen in 
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Fig. (J7|), the effective mass for SLy4 is larger than that of SkI3, we can understand why the 
entropy associated to SLy4 (left panel) is larger than the one corresponding to SkI3 (right 
panel) . 

Finally, using the low temperature expansion, we will analyze why pf decreases with 
temperature. At zero temperature, the difference between the free energy (internal energy) 
of the fully polarized phase and the non-polarized one is: 



A/(p,T = 0) = f(p,T = 0,A = l)-f(p,T 



5 2m*(p, A 
1 



-(3vr 2 p) 



2/3 



t (l-x ) + -t 3 (l 
o 



2 2/3 



0, A = 0) = 
m*(p, A 



m*(p, A = 0) 



(36) 



In order to find pf(T = 0), we should find the density at which the previous expres- 
sion becomes zero. For most of the forces, as we have already shown in Table |Hj the 
density at which the first term becomes negative is quite low, while the combination 
to(l — Xq) + 5*3(1 — x 3 )p a is usually negative at all densities. The energy difference, then, 
becomes negative only through a balance of both terms, and this happens at p = pf(T = 0), 
which is an upperbound for the real critical density p c (T = 0). It is quite noticeable that 
this cancellation results from a subtle interplay of all parameters of the Skyrme force. Once 
we introduce a low temperature, the difference between free energies becomes: 



f(p, T, A = 1) — f(p, T, A = 0) = A/(p, T = 0) 



7r 2 T2 2m*(p,A = 0) 
4 ft 2 (67r 2 p) 2/3 



m*(p, A 



2 2/3 



(37) 



m*(p, A = 0) 

and, therefore, this relation shows that at finite temperature, whenever m*(p, A = 
l)/m*(p, A = 0) > 2 2 / 3 , which occurs for densities larger than pf , the unpolarized system is 
unstable if the instability happens already at T = 0. However, even when A/(p, T = 0) > 0, 
the second term can make the difference become negative. Around pf(T = 0), the differ- 
ence decreases with density (see Fig. (|2*|) and thus the negative term makes pf (T 7^ 0) 
smaller with respect to pf (T = 0). Or, in other words, pf (T) < pf (T = 0) provided that 
Pc(T = 0) > pf, which is always true as seen in Table ITT1 



15 



V. SUMMARY AND CONCLUSIONS 



In this work, we have studied the properties of polarized neutron matter, with neutrons 
interacting through Skyrme-type interactions, both at zero and finite temperature. Firstly, 
we have revised the zero temperature calculations using two modern Skyrme forces and we 
have studied the ground state of neutron matter as a function of polarization. We have 
shown that the ground state is not necessarily at fully polarized or unpolarized matter, but 
that it can be found at partially polarized matter, giving rise to a spontaneous symmetry 
breaking where the system prefers a state with non-zero polarization. 

Our main emphasis, however, has been the study of the influence of temperature on the 
manifestation of the ferromagnetic behaviour. In particular, we have considered the stability 
of the unpolarized phase and we have shown that the critical density at which ferromagnetism 
takes place p c decreases with temperature. This unexpected behaviour has been associated 
to an anomaly of the entropy: above a certain critical density pf , the entropy of the polarized 
phase turns out to be larger than that of the unpolarized one. We have also shown that this 
fact is a consequence of the dependence of the entropy on the effective mass of the neutrons 
with different third spin component and, in particular, a consequence of the dependence 
of these effective masses on the polarization. More precisely, we have derived a condition 
for the maximum ratio between the effective masses of the fully polarized phase and the 
unpolarized one. Although this criterium is density dependent, one could in principle use it 
as a restriction on the parameters defining a Skyrme force. 

Finally, we would like to emphasize the fact that the present analysis has been restricted to 
Skyrme interactions which, contrary to other microscopic calculations, give a ferromagnetic 
transition at densities around 3.5po- It would be useful to compare these calculations with 
the results obtained with realistic interactions and determine the behaviour of the entropy 
and the effective masses of neutrons as a function of the spin polarization. Work in this 
direction is presently in progress. 
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Poo (fm" 3 ) 


a v (MeV) 


a s (MeV) 


^oo (MeV) 


M max (M ) 


SLy4 


0.160 


-15.97 


32.04 


230.9 


2.04 


SkI3 


0.158 


-15.98 


34.89 


259.2 


2.19 



TABLE I: Symmetric Nuclear Matter properties at saturation density for SLy4 and SkI3. The 
corresponding value of the maximum mass of a neutron star obtained with the EoS of these two 
forces is also shown. 
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FIG. 1: Ratio between the magnetic susceptibility of the free Fermi gas and the corresponding 
magnetic susceptibility of interacting neutron matter as a function of density for several tempera- 
tures. 
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pf (fm 3 ) 


p c (fm" 3 ) 


SGI 


0.09 


0.28 


SLyO 


0.11 


0.41 


SLyl 


0.15 


0.58 


SLylO 


0.13 


0.61 


SLy2 


0.09 


0.29 


SLy230a 


0.10 


0.54 


SLy3 


0.15 


0.60 


SLy4 


0.14 


0.60 


SLy5 


0.15 


0.57 


SLy6 


0.14 


0.90 


SLy7 


0.14 


0.57 


SLy8 


0.15 


0.60 


SLy9 


0.11 


0.47 


SV 


0.15 


0.77 


SkI3 


0.08 


0.37 


SkI5 


0.07 


0.28 



TABLE II: Densities pf at which the condition (|30|) is violated for 16 different Skyrme forces. The 
respective critical densities p c for the onset of ferromagnetism at T=0 are also shown. 



21 




FIG. 2: Difference between the free energy per particle of fully polarized neutron matter and 
unpolarized neutron matter as a function of density for several temperatures. 
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FIG. 3: Difference between the energy per particle of fully polarized neutron matter and unpolarized 
neutron matter as a function of density for several temperatures. 
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FIG. 4: Difference between the entropy per particle of fully polarized neutron matter and unpo- 
larized neutron matter as a function of density for several temperatures. 
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FIG. 5: Left panel: neutron matter free energy per particle at zero temperature as a function of 
polarization for several densities. Right panel: neutron matter free energy per particle at a fixed 
density p = 0.36 fm -3 as a function of polarization for several temperatures. Both figures were 
obtained using the SkI3 force. 
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FIG. 6: Polarization associated to the minimum of the free energy for neutron matter as a function 
of density for several temperatures. 
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FIG. 7: Neutron effective mass of the up and down components as a function of polarization at 
p = 0.16 fm -3 . The horizontal lines signal the maximum value of the effective mass of the up 
component in fully polarized matter if the entropy of the polarized phase has to be smaller than 
the unpolarized one. 
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FIG. 8: Entropy per particle of neutron matter as a function of polarization at p = 0.32 fm 3 for 
several temperatures. 
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